Static disorder and local structure in zinc(II) isonicotinate, 
a quartzlike metal— organic framework 



Ines E. Collings, 1 Matthew G. Tucker, 2 David A. Keen, 2 and Andrew L. Goodwin 1 ^ 

1 Department of Chemistry, University of Oxford, 
Inorganic Chemistry Laboratory, South Parks Road, Oxford 0X1 3QR, U.K. 
2 ISIS Facility, Rutherford Appleton Laboratory, Harwell Science and Innovation Campus, Didcot 0X11 OQX, U.K. 

(Dated: December 30, 2011) 

Using a combination of Rietveld and RMC refinement of neutron total scattering data, we find that 
the 10 K structure of zinc(II) isonicotinate shows strong evidence of static disorder. This disorder 
takes the form of transverse displacements of the isonicotinate ligand and results in elongated atomic 
displacement parameters and dampened oscillations of the experimental G(r). We analyse the 
RMC configurations using an approach derived from geometric algebra. Complications regarding 
the inclusion of hydrogenous guest molecules within the pore structure are discussed. This study 
highlights the way in which structural flexibility can give rise to multiple low-energy ground states 
in MOF-type materials. 



I. INTRODUCTION 

The structural analogies between traditional oxide- 
containing frameworks and the class of materials known 
as metal-organic frameworks (MOFs) are certainly well 
established. 1 A topical example is the family of zeolitic 
imidazolate frameworks (ZIFs), which are MOF ana- 
logues of zeolitic Si02 polymorphs. 2 " 4 Replacing Si 4+ and 
2 ~ ions by tetrahedral Zn 2+ centres and bridging imi- 
dazolate linkers respectively, it is possible to generate a 
wide range of zeolite-like framework structures with in- 
creased pore volume. Like zeolites, ZIFs offer attractive 
gas storage and catalytic functionality,^ but they couple 
these properties with the characteristic chemical versa- 
tility for which MOFs are renowned (e.g. the imidazolate 
linker can be replaced by substituted derivatives and the 
transition metal dication varied) . The extensive range of 
functionalities exhibited by ZIFs and many other MOF 
families has generated a strong and sustained interest 
in studying structure/property relationships in such 'hy- 
brid' inorganic-organic materials}^! consequently, there 
is a great deal now known regarding their structural 
chemistry. 

What has become increasingly evident is that the 
structural behaviour of MOFs is often significantly more 
complex than might have been anticipated. A sur- 
prisingly large number — including the canonical MOF- 
5 (Ref. [TUj) — exhibit negative thermal e xpans ion (NTE; 
i.e., their lattices expand on cooling) I 11 ! 12 ! Some un- 
dergo pressure- and temperature-induced amorphisation 
processes.^ Yet others show unprecedented structural 
flexibility, being capable of dramatic changes in crystal- 
lite dimensions during absorption or desorption of guest 
species.^ While it is the case that similar mechanical 
phenomena have been documented for oxide-containing 
frameworks (e.g. NTE in ZrW 2 8 , see Ref. US), the ef- 
fects observed in MOFs are almost always much more 
extreme. The increased structural flexibility of molecu- 
lar metal-ligand-metal linkages is heavily implicated in 
the fundamental difference in magnitude of behaviour: 



metal-ligand geometries are more easily distorted when 
the ligand is large and flexible than when it is a single 
oxide ion. 

The prevalence of low-energy deformation mechanisms 
in MOFs should result in an increased propensity for 
structural distortion, which if incoherent could resem- 
ble the static disorder often found in zeolites^ and per- 
ovskite frameworks^ — or when coherent may give rise to 
an extreme symmetry lowering observed for some com- 
plex framework oxidesP^ Yet it seems that reports of 
static disorder in MOFs are relatively few in number: 
there is, for example, some discussion of ligand orienta- 
tion disorder and strongly anisotropic displacement pa- 
rameters in Cu(4-oxopyrimidinate)2p^ likewise molecu- 
lar dynamics studies of MOF-5 and its isoreticular con- 
geners have also pointed to the existence of low-barrier 
enthalpy landscapes though t to b e responsible for both 
static and dynamic disorder! 20 ! 21 ! A traditional emphasis 
within the field on average structure studies — which are 
also usually performed at temperatures insufficiently low 
to distinguish static disorder from vibrational motion — 
has perhaps meant that the degree of static disorder 
within MOFs has remained poorly understood. 

Here, as part of a broader study into the local struc- 
ture and dynamics of mincralomimetic MOFs, we re- 
port a combined average- and local-structure investiga- 
tion into the existence and nature of static disorder in 
zinc(II) isonicotinate, Zn(ISN)2- Our approach has been 
to collect neutron total scattering data for Zn(ISN)2 at 
10 K, which we analyse using Rietveld refinement (av- 
erage structure) and reverse Monte Carlo (RMC) re- 
finement of the corresponding pair distribution function 
(PDF) data (local structure). We find strong evidence in 
both types of refinement for the existence of static disor- 
der that resembles local transverse displacements of the 
isonicotinate ligands. 

The Zn(ISN)2 framework structure contains Zn 2+ cen- 
tres that are fourfold-coordinated by isonicotinate an- 
ions, with each anion bridging two Zn 2+ cations to form 
a tetrahedral net with the quartz topology [Fig. [I] . The 
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FIG. 1: Crystal structure of (a) Zn(ISN) 2 with [ZnN 2 2 ] 
tetrahedra connected via isonicotinate ligands, and (b) a- 
quartz with connected [SiCU] tetrahedra. Both structures are 
viewed down the c axis with their unit cell shown. Local 
coordination environment of (c) Zn(ISN) 2 and (d) a-quartz. 
(Zn, green; N, light blue; O, red; C, black; Si, dark blue, and 
isonicotinate hydrogens have been removed for clarity). 



coordination sphere of each Zn 2+ centre consists of two 
pyridinyl N donor atoms and two carboxylate O atoms — 
an arrangement that is inconsistent with the correspond- 
ing site symmetry imposed by the /3-quartz structure. In- 
stead, the material adopts the lower-symmetry a-quartz 
structure, crystallising in the hexagonal space group 
P6 2 !22l As in quartz itself, the structure of Zn(ISN) 2 
contains one-dimensional pores parallel to the crystal- 
lographic c axis. In as-prepared samples these pores are 
occupied by solvent molecules, but the pores can be evac- 
uated by heating to 100 °CP There remains some contro- 
versy as to whether the correct space group for Zn(ISN)2 
is actually P3i , which is a subgroup of -P62 and which de- 
mands the existence of two cryst allog raphically distinct 
ISN units in the asymmetric unit 122123] 



Our paper begins by describing the experimental meth- 
ods used in our study, together with our approach to 
RMC refinements of the neutron PDF data we have col- 
lected. We note that this is somewhat of a test-case for 
the use of our RMCPROFILE code^ in the study of struc- 
tural disorder in crystalline MOFs, and we have encoun- 
tered some difficulties through an unquantifiable degree 
of absorption of atmospheric H2O into our deuterated 
sample. Mindful of the ensuing limitations, we present 
the results of our average-structure and RMC studies and 
discuss the influence of the included solvent on the con- 
clusions drawn. We finish with a discussion on the nature 
of static disorder in this material in the context of corre- 
lated (thermal) displacements within cv-quartz itself. 



II. MATERIALS AND METHODS 

A. Synthesis 

A polycrystalline sample of zinc(II) isonicotinate was 
prepared by mixing stoichiometric quantities of zinc (II) 
acetate and d 4 -isonicotinic acid (QMx, 98% D) dis- 
solved in dimethylsulfoxide (DMSO). The white powder 
formed on mixing was filtered, washed and dried in vacuo 
(100 °C, 24 h) in order to remove all solvent from the 
framework pores. Working within a glovebox, the dried 
sample (ca 2g) was transferred to a vanadium can suit- 
able for neutron total scattering measurements. 



B. Neutron total scattering 

Neutron total scattering data were collected at 10 K 
for Zn(ISN) 2 using the time-of-flight diffractometer GEM 
at ISISpSHlZl p or the experiment, approximately 2g 
of Zn(ISN) 2 , prepared as described above, was placed 
within a cylindrical vanadium can of 8.3 mm diameter 
and 6 cm height. This can was loaded at room temper- 
ature inside a closed cycle helium refrigerator and the 
temperature was lowered slowly to 10 K. Total scattering 
data were collected over a large range of scattering vec- 
tors of magnitudes 0.6 < Q < 40 A -1 , giving a real-space 
resolution of order Ar ~ 0.09 A. 

The total scattering data were corrected using stan- 
dard methods, taking into account the effects of back- 
ground scattering, absorption, multiple scattering within 
the sample, beam intensity variations, the Placzek inelas- 
tic correction, and hydrogen content corrections.^ These 
corrected data were then converted to experimental G(r) 
and F(Q) functionsPHH 



F(Q) = po I 4^r 2 G(r)^dr 
la Q r 



G(r) 



CiCjbibjigijir) - 1], 
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riij(r) is the number of pairs of atoms of type i and j 
separated by distance r, po is the number density, Cj the 
concentration of each species i and hi the corresponding 
neutron scattering length. The Bragg profiles for each 
data set were extracted from the scattering data collected 
using the detector banks centred on on scattering angles 
26» = 34.96°, 63.62°, 91.30°, and 154.40°. 

Because the number density p enters Eqs. ([!]) and 
it is possible during data normalisation to check for con- 
sistency between the expected value of po and the value 
for which normalisation is most robust. In the present 
case, a value of po corresponding to guest-free Zn(ISN) 2 
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was found not to be the value most consistent with the 
observed data; moreover, by comparing the incoherent 
scattering levels associated with each detector bank it 
was possible to deduce that the sample as measured was 
actually contaminated with a significant quantity of hy- 
drogenous material, despite the care taken to evacuate 
the framework completely. 

This unquantifiable degree of solvation had three con- 
sequences for our (usually quantitative) RMC modelling 
of the PDF data. First, we found that quantitative nor- 
malisation of the G(r) function was not possible, and 
so for RMC refinement a smoothly-varying function was 
subtracted from the normalised data — this did not af- 
fect the final model but improved visually the quality of 
our fits to data. Second, it was not possible to fit the 
very lowest-r region of the PDF since this was likely to 
contain a well-structured contribution from the included 
solvent (e.g. the O-H distance if the solvent were H2O). 
Third, because the value of po was n °t known accurately 
it was necessary to fit the PDF only while allowing re- 
finement of an overall scale parameter. We note that 
such an ap pro ach is quite common for other PDF fitting 
procedures 



C. Average structure refinement 

The experimental Bragg diffraction profiles wer e fit- 
ted with the GSAS Rietveld refinement progranj 30 * 31 * using 
both the published P3i and P6 2 structural modelsp2E31 
Atomic coordinates of the organic ligand were refined 
using rigid body constraints in order to minimise the 
number of refinable parameters; likewise a single set of 
anisotropic displacement parameters was refined for all 
atoms in the isonicotinatc group. Zn atom coordinates 
were refined while fixing the z component (there is no 
unique origin in z for either P3i or P62) and an isotropic 
displacement parameter was allowed to refine freely. Sol- 
vent occupancy within the pore network was treated us- 
ing a number of different approaches, and these are dis- 
cussed in more detail in the results section below. 



D. Reverse Monte Carlo refinement 

RMC refinements were carried out using the rmcpro- 
file codePS To the best of our knowledge, this study 
represents the first RMC study of a crystalline MOF and 
we detail below some of the difficulties we have encoun- 
tered in the process. As for all RMC studies, the basic re- 
finement objective is to produce atomistic configurations 
that fit simultaneously the experimental G(r), F(Q), and 
Bragg profile I(t) functions. This is achieved by accept- 
ing or rejecting random atomic moves produced by the 
Metropolis Monte Carlo algorithm, where in this case the 
Monte Carlo 'energy' function is determined by the qual- 
ity of the fits to data. The motivation behind fitting both 
real-space G(r) and reciprocal-space F(Q)/I(t) functions 



TABLE I: 'Distance window' parameters d m i n , d max used in 
our Zn(ISN)2 RMC refinements and the corresponding mean 
pair separations d (A) and their standard deviations a. 



Atom pair 


t^min (A) 


G^max (A) 


J (A) 


a (A) 


O-H 


0.8 


l.l 


1.006 


0.076 


C-D 


0.9 


1.2 


1.039 


0.072 


C-C 


1.1 


1.7 


1.431 


0.062 


C-N 


1.1 


1.7 


1.321 


0.055 


C-0 


1.1 


1.7 


1.268 


0.040 


Zn-0 


1.8 


2.4 


2.004 


0.128 


Zn-N 


1.8 


2.4 


2.081 


0.104 



is to probe local distortions in the framework in a manner 
that is inherently consistent with the long-range periodic 
order reflected in the Bragg intensities. Similar refine- 
ments of crystalline materials in which reciprocal-space 
data are not used have a tendency to become unphysi- 
cally disordered,^ and we were keen to assess the degree 
of structural disorder in Zn(ISN) 2 on the most realistic 
level possible. 

Our starting configurations for the RMC process were 
based on a orthogonal supercell related to the crystallo- 
graphic cell by the transformation 



a 




"4 " 




a 


h 




2 4 


X 


b 


c 


RMC 


10 




c 



P62 



giving a cell of dimensions 60.14 A x 53.82 A x 61.36 A. 
The use in RMC of orthogonal axes for hexagonal sys- 
tems means the configurations can be prepared with di- 
mensions approximately equal in each direction. This 
maximises the pair distribution cut-off value r max for a 
given number of atoms (and hence minimises computa- 
tional cost). 

In the results section below we discuss in more detail 
the various tests we performed to determine the best way 
of modelling solvent occupancy. The key set of RMC 
configurations contained 18 720 atoms, of which 12 960 
are part of the framework and 5 760 are part of the in- 
cluded solvent. We note that the configuration contained 
six different atom types (C, D, H, N, O, Zn), giving rise 
to a total of 21 partial gij(r) contributions — we under- 
stand this to be the largest number refined within RM- 
CPROFILE to date. A number of soft constraints and re- 
straints were applied throughout the refinement process: 
closest-approach constraints stopped atom pairs from be- 
ing separated by unphysically short distances; 'distance- 
window' constraints maintained the framework connec- 
tivity without prejudicing the bond-bending and bond- 
stretching terms [Table[l]p^and, finally, a set of empirical 
geometric restraints were applied to maintain the geom- 
etry of the isonicotinate moietyP^ 

The refinement process included fitting the experimen- 
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tal F(Q) (0.6 < Q < 40 A" 1 ), I(t) and G(r) (1.6 < r < 
25.2 A) functions. We note here that the minimum value 
of r used for the G(r) fitting is larger than the nearest- 
neighbour bond length for the vast majority of solvents, 
including H2O. A total of five independent RMC refine- 
ments were performed in parallel; each refinement was 
allowed to continue until no further improvements in the 
fits to data were observed. The absolute atomic coor- 
dinates differ amongst the five final configurations, but 
the corresponding fits to data are essentially identical. 
Wherever possible, our results are averaged over all five 
configurations. 



III. RESULTS 
A. Average structure 

We used as the basis of our Rietveld refinements the 
models of Refs. [2U and [231 both of which include a num- 
ber of water molecules within the framework pore struc- 
ture. Our first set of refinements involved simply remov- 
ing the solvent component and refining framework coor- 
dinates. We obtained essentially identical fits for both 
P3i and P62 models, neither of which was entirely con- 
vincing [Fig. [2]Ja)]. A considerable improvement to the 
quality of these fits was obtained by including scattering 
density within the framework pores. We tested a number 
of different ways of modelling this scattering density — 
including the coordinates of Refs. |2"21 and |2"31 — but found 
the best fits were obtained for a slightly different struc- 
ture model containing four H2O molecules per formula 
unit [Fig. |2jb)]. The corresponding crystallographic de- 
tails are summarised in Table [TTJ 

It is certainly feasible that the sample absorbed atmo- 
spheric H 2 during setup of the GEM experiment or that 
some H2O or DMSO (used in the synthesis) remained 
from incomplete desolvation during sample preparation. 
Because all of the various models for adsorbed solvent we 
tested involved large displacement parameters and frac- 
tional occupancies of multiple equivalent sites, we are 
certainly not claiming that the diffraction intensities are 
sufficiently sensitive to identify the chemical composition 
of the included component. Instead, we were keen to 
assess the extent to which the structural parameters of 
the framework itself were affected by the different mod- 
els. The positions and atomic displacement parameters 
obtained for the 'four H2O' and solvent-free models are 
given in Table [TTJ their close correspondence illustrates 
that the framework geometry is robustly defined by the 
diffraction data, even if those of the solvent component 
are not. In all our refinements, we found no evidence to 
support symmetry lowering to the P3i model of Ref. 1231 

It was possible to refine reliably a set of anisotropic dis- 
placement parameters for the ISN ligand. The values ob- 
tained correspond to thermal ellipsoids that are strongly 
elongated in a direction perpendicular to the plane of 
the ligand itself [Fig. |3](a)]. This sort of behaviour is not 
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FIG. 2: Representative Reitveld fits to the 10 K neutron 
diffraction pattern of Zn(ISN)2: (a) using an evacuated frame- 
work model, and (b) using a model containing four water 
molecules per formula unit. Experimental data are given as 
small filled circles, the fitted profile is shown as a solid red line, 
and the difference (fit— data) is shown beneath each curve. 

uncommon for molecular framework materials, since the 
lowest-energy vibrational modes usually involve trans- 
verse displacements of the bridging ligand; examples in- 
clude the cyanide bridges of Zn(CN) 2 (Ref. [36]) and the 
terephthalate linkers in MOF-5 fRef. Wf). What is un- 
usual here is the magnitude of the thermal ellipsoids 
given that the data were collected at 10 K. We note that 
by using a wide range of GEM detector banks our refine- 
ments include data at sufficiently low d-spacing values to 
give confidence in the determination of anisotropic dis- 
placement parameters. The issue of the magnitude of 
displacement parameters is one to which we return in 
Section El 

B. Local structure 

With the strong preference in Rietveld refinements for 
a structure model containing four water molecules per 
unit cell, our RMC study of local structure in Zn(ISN) 2 
used this same model as its starting point. RMC refine- 
ment gave fits to the time-of-flight Bragg intensity func- 
tion I(t) [Fig. [4] that were essentially identical to those 
obtained by Rietveld refinement. As anticipated, the sol- 
vent molecules were strongly disordered in our final RMC 
configurations and consequently our analysis focusses on 
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TABLE II: Crystallographic parameters, atomic coordinates 
and isotropic equivalent displacement parameters determined 
using Rietveld refinement of neutron scattering data for 
Zn(ISN)2. The ISN coordinates are given as the centre-of- 
mass of the ligand. The first set of positional and displace- 
ment parameters correspond to the final structural model 
containing four H2O molecules per formula unit; the sec- 
ond (marked with an asterisk) correspond to the solvent-free 
model discussed in the text. 



Crystal system Hexagonal 



Space group P62 
a (A) 15.53615(25) 
c (A) 6.13556(26) 
V (A 3 ) 1282.54(6) 
Z 3 

T (K) 10 

Atom x y z C/iso (A 2 ) 

Zn 0.5 0.1319(16) 0.0199(22) 

ISN 0.26272(8) 0.79764(9) 0.4963(12) 0.0229(3) 

H 2 01 0.7702(11) 0.8083(6) 0.310(4) 0.438(5) 

H 2 Q2 0.8143(8) 1.0246(7) 0.3744(20) 0.438(5) 

Zn* 0.5 0.132(3) 0.016(4) 

ISN* 0.26305(15) 0.79501(17) 0.514(3) 0.0204(7) 



a 




FIG. 3: Unit cell projections of Zn(ISN) 2 viewed down the 
c axis, showing (a) the average structure and (b) the atomic 
distributions obtained by collapsing the RMC configuration 
onto a single unit cell. In both cases, water molecules have 
been removed for clarity. 



the local structure of the Zn(ISN) 2 framework itself. 

By 'collapsing' the atomic coordinates of each refined 
RMC configuration onto a single unit cell, it is possible 
to visualise the average structure model to which these 
configurations correspond [Fig. [3j»]. There is a clear 
similarity to the results of our Rietveld refinement. In 
particular, there are large anisotropic displacements of 
the ISN ligand that reflect well the displacement param- 
eters determined above. Whereas the Rietveld refine- 
ment employed rigid-body constraints and a single set 
of displacement parameters for all atoms within the ISN 
ligand, the RMC refinement allows the scattering dis- 
tribution for each atomic site to assume whatever form 
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FIG. 4: RMC fit to the Bragg profile with the experimental 
data points given as points, and the RMC fits (solid lines) 
obtained using RMCPROFILE as described in the text. The 
difference curve is shown beneath the fit. 
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FIG. 5: Calculated and experimental G(r) functions for 
Zn(ISN) 2 . (Bottom) Experimental G(r) data (filled black 
circles) and RMC fit (red line). (Middle) Calculated G(r) 
function for the static disordered model obtained from Ri- 
etveld refinement. (Top) Calculated G(r) function for the 
low-energy-dynamics model described in the text. 

is demanded by the data. We find that the differences 
amongst distributions for the various atoms in the ISN 
ligand are small, and certainly in all cases the displace- 
ments are much stronger in a direction perpendicular to 
the ligand plane. 

Given the ambiguity regarding included solvent com- 
position, it is not at all surprising that we found the 
lowest-r region of the G(r) particularly difficult to nor- 
malise robustly. For most of the G(r), the contribution 
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from a disordered solvent network would be expected to 
be a smoothly- varying function of r; however, at the low- 
est values of r there must exist well-defined contributions 
corresponding to the interatomic separations within indi- 
vidual solvent molecules. Consequently our RMC fits to 
G(r) were performed only for r > 1.7 A in order that the 
particular method employed for modelling solvent inclu- 
sion (in this case four H2O molecules per unit cell) would 
not affect refinement of the framework. The fits obtained 
are not at the quantitative level usually achieved using 
rmcprofile, but nonetheless all qualitative features of 
the G(r) are certainly captured well [Fig. [H]. 

We were struck by the absence of pronounced features 
in the PDF beyond r ~ 5 A, and especially so for data 
collected at a temperature of 10 K. Even for a system with 
a disordered component, the framework itself would be 
expected to give rise to a well-structured G(r) function: 
the superposition of two functions — one strongly vary- 
ing and the other weakly varying — still varies strongly. 
Consequently, the G(r) function could be considered con- 
sistent with a large degree of disorder in the framework 
geometry of Zn(ISN)2. 

It is straightforward to determine bond-angle distri- 
butions directly from the RMC configurations, and we 
concentrate here on two particular distribution functions. 
The first corresponds to the angles within the [ZnN202] 
coordination tetrahedra [Fig. p[a)] . The second type of 
angle is based on the centres of mass of the ISN ligands 
(which we represent by the symbol 'X'): then the X-Zn-X 
angle distribution [Fig. |6jb)] is analogous to the O-Si- 
O tetrahedral angle distribution in a-quartz itself. We 
find in both cases that the distributions are very broad, 
suggesting that Zn(ISN)2 is flexible both on the scale of 
the framework itself and also in terms of the individual 
coordination polyhedra. 

The degree of structural distortion can be quantified 
further using the language of geometric algebra (GA)P^ 
Analysis of RMC configurations using the GA-based code 
GASpMfini allows the atomic displacements to be under- 
stood in terms of (i) translations, (ii) rotations, and (iii) 
distortions of fundamental geometric units. In this par- 
ticular system, there are two types of geometric unit of 
special interest. The first corresponds to the [ZnN202] 
coordination tetrahedra, the behaviour of which reflects 
the rigidity of Zn-N and Zn-0 bonding interactions. The 
second type of unit involves the centres of mass of the 
ISN ligands: here the behaviour of the [Z11X4] tetrahedra 
describes the flexing and deformation of the framework 
structure as a whole. 

In Fig. [7] we plot the relative components of transla- 
tions, rotations, and deformation (stretching and bend- 
ing) determined for the two types of polyhedral unit 
across our five RMC configurations. What we find is 
that the translational component is very large indeed for 
the Zn coordination polyhedra. These translations are 
evidently allowed by flexing of the Zn-ISN-Zn linkages, 
reflected in the large angle bending component of the 
[ZnX4] values. These distortions will correspond to in- 
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FIG. 6: RMC bond angle distributions: (a) intratetrahe- 
dral angles within the [ZnN2 02] coordination environment 
(N-Zn-N in blue, O-Zn-0 in red, and N-Zn-0 in black); (b) 
intratetrahedral X-Zn-X angles reflecting geometric flexing 
of the framework itself. 
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FIG. 7: Average mismatch scores obtained using GASP for 
a-quartz (black), the [Zn02N2] coordination environment in 
Zn(ISN)2 (dark grey), and the [ZnX.4] framework geometry 
of Zn(ISN)2 (light grey). The bar chart shows the type of 
distortions present in the polyhedra: (left-right) translations, 
rotations, bond angle bending and bond stretching. 



creased displacement of the Zn atoms away from their 
average sites, and also an apparent increase in the trans- 
verse component of the anisotropic displacement param- 
eters for the ISN ligands. These two effects are precisely 
those observed in Rietveld refinement. 

Finally, we put these results in context by comparing 
the magnitude of the translation, rotation and distor- 
tion components with those obtained for an RMC refine- 
ment of a-quartz itself (also at 10 K) [Fig. [7] Not only 
are the rigid-body-type translations and rotations many 
times larger for Zn(ISN)2, but the distortions are also an 
order of magnitude more extreme. This result is consis- 
tent with the general perception that MOF-type materi- 
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als are substantially more flexible than their oxide-based 
counterparts. 



IV. DISCUSSION 

That the ISN ligands in Zn(ISN)a are displaced in 
a transverse direction to a very large degree even at 
10 K seems clear both from the average- and local- 
structure analysis performed so far. The central ques- 
tion is whether or not the magnitude of this displace- 
ment is consistent with low-energy vibrational motion of 
the linkages or with static disorder. Here we make use 
of the following lattice dynamical formalism which links 
magnitude of displacement with phonon mode energies: 



- +n(w E ,T) 



(5) 



where n(o>E,T) = 1 / {exp[hu)E / feT] — 1} is the Bose- 
Einstein occupation number and m the atomic mass. 

Because 10 K is low with respect to typical phonon fre- 
quencies, the displacements in equation Eq. ^ will be 
dominated at this temperature by the energy of the first 
dispersionless branch of the phonon spectrum. Phonon 
frequencies have not been determined for Zn(ISN)2, but 
we might expect that the relevant span of energies will 
be roughly similar to that in other MOFs; e.g. MOF- 
5, the structure of which is also based on Zn 2+ cen- 
tres but connected via terephthalate ligands (similar in 
size to isonicotinate) The first dispersionless branch 
in MOF-5 is calculated to occur at uj = 4.6 THz. Using 
this value as a conservative estimate for we and substi- 



tuting into Eq. (|5fo one obtains 



4 Zn/ 



0.0113 A 2 and 
( m isn) — 0.0059 A 2 . Comparing these values to the Ui so 
values in Table [TTJ it is clear that the experimental Zn 
and ISN displacements are, respectively, two and four 
times larger than can be accounted for by thermal mo- 



tion alone. This is clear evidence that the displacements 
observed for Zn(ISN)2 are consistent only with static dis- 
order of the framework. 

As a further check of this conclusion, we calculated us- 
ing PDFGUI (Ref. |43|) the G(r) functions expected from 
the Rietveld model of Section |III A| using both the as- 
refined Uij/Ui so parameters (i.e. with static disorder) 
and using equivalent parameters re-scaled according to 
the (uf) values estimated above (i.e. for a comparable 
model where displacements are due only to low-energy 
vibrational modes) . The results of both calculations are 
shown in Fig. [5j from which it is clear that (i) the dy- 
namic model gives a PDF that is much more typical of 
10 K data, and (ii) the static disorder model reflects well 
the actual experimental G(r). 

In conclusion, we have used a combination of Ri- 
etveld and RMC refinement of neutron total scattering 
data to study the local and average structure of the 
metal-organic framework zinc(II) isonicotinate. While 
our study has unquestionably been complicated by the in- 
clusion of hydrogen-containing solvent molecules within 
the pore structure of the framework, we do find robust 
evidence for the existence of large-scale static disorder 
in the material at 10 K. The role of structural flexibility 
in the driving unusual dynamic properties of MOFs is 
increasingly appreciated; what we find here is that the 
same flexibility may also affect the ground-state struc- 
tural properties of the same materials. 
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